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^— ^ ' Abstract: We study the dynamics of planetary systems with two planets moving in the same plane, 

when frictional forces act on the two planets, in addition to the gravitational forces. The model of the 

C^ I general three-body problem is used. Different laws of friction are considered. The topology of the phase 

space is essential in understanding the evolution of the system. The topology is determined by the families 

of stable and unstable periodic orbits, both symmetric and non symmetric. It is along the stable families, 

or close to them, that the planets migrate when dissipative forces act. At the critical points where the 

stability along the family changes, there is a bifurcation of a new family of stable periodic orbits and 

n , ' the migration process changes route and follows the new stable family up to large eccentricities or to a 

["tI I chaotic region. We consider both resonant and non resonant planetary systems. The 2/1, 3/1 and 3/2 

(-H ' resonances are studied. The migration to larger or smaller eccentricities depends on the particular law 

f^ I' of friction. Also, in some cases the semimajor axes increase and in other cases they are stabilized. For 

Q , particular laws of friction and for special values of the parameters of the frictional forces, it is possible 

to have partially stationary solutions, where the eccentricities and the semimajor axes are fixed. 

C/2 '. 

C^ . Keywords: Periodic orbits, resonances, extrasolar planetary systems, migration. 
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.^ ! 1 Introduction 

in 

^^ I Since the last decade of the 20th century it is known that there are 340 observed extrasolar 

• ' planetary systems with more than 400 planets and some of them have two or more planets 

^— ^ . (Schneider, October 2009). In many cases of multiplanetary systems, two planets are in mean 

f^ I motion resonance, (e.g. HD 82943, GJ 876, HD128311, 55CnCb,c)- Some of these systems have 

large eccentricities and are evidently stable, since they exist in nature. Thus the role of the 
resonances in the evolution of a planetary system should play an important role either for its 

k> i long term stability or its migration to the present position. 

j^ ■ There are different approaches to the study of the long term dynamical stability of a resonant 

planetary system and on the mechanisms that stabilize the system or generate chaotic motion 
and instability (e.g. Beauge et al, 2003; Ferraz-Mello et al, 2006; Marzari et al, 2006; Callegari 
et al, 2006; Baluev, 2008; Gozdziewski et al, 2008; Beauge et al, 2008; Schward et al, 2009). 
In these papers different methods have been applied, as the averaging method, direct numerical 
integrations of orbits, or various numerical methods which provide indicators for the exponential 
growth of nearby orbits. In this way the regions where stable motion exists have been detected, 
in the orbital elements space. Another method to study the dynamics of a planetary system 
is to find the basic families of periodic orbits and their stability. It is known that the periodic 
orbits play a dominant role in understanding the dynamics of a system, because they determine 
critically the structure of the phase space. In this way, we can detect the regions where stable 
librations could exist (Hadjidemetriou, 2006; Voyatzis, 2008). These are the regions close to 
the stable periodic orbits, where a real planetary system could be trapped. The periodic orbits 
with nonzero eccentricities of the planets, either stable or unstable, are associated with mean 
motion resonances between the planets. Also, work on these lines are by Hadjidemetriou (2002), 
Psychoyos and Hadjidemetriou (2005), Voyatzis and Hadjidemetriou (2005, 2006). A review on 
these issues is given in Hadjidemetriou (2008). 



It is believed that the existence of a planetary system in a region of long term stability is the 
result of a dynamical process that forces the planets to migrate after their formation (see e.g. 
Tsiganis et al, 2005; Morbidelli et al, 2007). Friction, due to the interaction between the newly 
formed planets and the protoplanetary nebula, before this latter is dissolved, is the main source 
of non conservative forces that generate a migration of the planets. Also tidal friction due to 
the body deformation also produces non conservative forces that affects both the orbit and the 
spin evolution (Ferraz-Mello 2008). Non conservative forces may explain the large eccentricities 
that are observed in some planetary systems or the very small planetary distances from the sun. 

The introduction of non conservative forces destroy the symplectic structure of a purely gravita- 
tional system and the phase space volume is not conserved. Consequently the overall dynamics 
and the topology of the phase space change critically. A particular review and relative literature 
can be found in Contopoulos (2002). 

Numerical simulations using hydro dynamical models (e.g. Nelson and Papaloizou, 2003a, b) or 
dissipative restricted three body models (Beauge and Ferraz-Mello, 1993; Gomes, 1995; Haghigh- 
ipour, 1999) showed that planetary evolution leads to a trapping of planets in a mean motion 
resonance. Such evolution follows paths of resonant stable configurations (Lee and Peal 2002; 
Papaloizou 2003; Lee, 2004; Sandor and Kley, 2006). Such paths are proved to be character- 
istic curves of resonant corotations (Ferraz-Mello et al, 2003; Beauge et al 2006; Zhou et al, 
2008), which actually coinside with families of periodic orbits in a rotating frame (Voyatzis and 
Hadjidemetriou 2005, 2006). The present work aims to study the close association of a possible 
evolution of a planetary system under dissipative forces with the families of periodic orbits of 
the conservative model. As we shall see in the following, the families of periodic orbits play a 
crucial role on the evolution of the system, in phase space, under any dissipation law and the 
knowledge of the families of periodic orbits is necessary in understanding the migration process. 

We consider, additionally to the gravitational forces, the action of the Stokes non conservative 
force 

R = -10-''{v-Vc), iy>0 (1) 

where the velocity v'c is the Keplerian circular velocity at the distance r from the Sun, given by 
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The vector eg is the unit vector normal to the radius vector r and rriQ is the mass of the sun. 
The Stokes dissipation law ([TJ implies that the nebula that generates the dissipation rotates 
differentially with a Keplerian circular velocity Vc at any distance r. We also consider a variation 
of the Stokes force given by 

R = -10-''{v-vJr^/^), iy>0. (2) 

This law differs from the law given in Equation ([T]), in that the circular velocity is multiplied 
by 1/r^''^. This means that the velocity of rotation of the protoplanetary nebula varies with the 
distance from the sun. It is larger than Vc close to the sun and goes slowly to zero as r — > cxd. 
We remark that R is applied to both planets and not only to the outer planet as e.g. in the 
studies of Lee and Peal (2002), Beauge et al (2006). 

The model we use is the general three body problem, for planar motion, with the sun, S, and the 
two planets. Pi (inner) and P2 (outer), as the three bodies of mass tuq, mi and m2, respectively, 
normalized such that mQ + mi + m2 = 1 . We assume that the center of mass of the whole system 
is at rest with respect to an inertial frame XOY. We have four degrees of freedom, for planar 
motion, with generalized variables Xi, Yi, X2, Y2. We define next, a rotating frame of reference 



xOy , whose x-axis is the hne S — Pi, with origin at the center of mass of these two bodies, where 
S is the Sun and Pi the inner planet. This is a non uniformly rotating frame, where the planet 
Pi moves on the x-axis and the planet P2 moves in the xy plane. As generalized coordinates 
in this rotating frame we can select xi, X2, y2 and 9, where 9 is the angle between the Ox axis 
and a fixed direction in the inertial frame. It turns out that the angle 9 is ignorable, so we have 
three degrees of freedom and the study of the system is restricted in the rotating frame xOy, 
with variables xi, X2-, 2/2 (Hadjidemetriou, 1975). 

The evolution of the planetary system can be studied by computing the Poincare map on a 
surface of section. By this method we reduce the dimensions of the phase space, without losing 
the generality of the problem. In the present study, we consider the surface of section 

y2 = (2/2 > 0), 

in the rotating frame xOy. 

It can be shown (see e.g. Voyatzis and Hadjidemetriou 2005, 2006) that families of symmetric or 
asymmetric periodic orbits exist in the rotating frame xOy. In particular, in a symmetric periodic 
orbit the planet P2 intersects perpendicularly the x-axis while the planet Pi is temporarily at 
rest on the x-axis. Consequently, the non zero initial conditions of a symmetric periodic orbit are 
2^10) 2:20) ^20) and a family of symmetric periodic orbits is represented by a smooth curve in the 
space 113 = {xio, X205 2/20} of initial conditions. An asymmetric periodic orbit is represented by a 
smooth curve in a five dimensional space of initial conditions e.g. lis = {xiQ,X2Q-,xiQ,X2Q-,y2o}- 
In the following we shall present a family by its projection on the coordinate plane X10X20 or in 
the space 6162 of the osculating eccentricities at t = 0. When we plot symmetric periodic orbits 
in the space 6162 we use the convention e^ > or Cj < when the planet is at aphelion or at 
perihelion, respectively. 

In the conservative system the stable fixed points of the Poincare map correspond to periodic 
orbits in the rotating frame. If in the equations of motion we consider dissipative forces, in 
addition to the gravitational forces, the energy integral does not exist. The stable fixed points 
are now stable limit cycles, to which the system is attracted in some cases. Chaotic attractors 
may also exist and correspond to a set of points (as t — > 00), which are irregularly distributed 
in the section. 



2 The evolution of a non resonant orbit under dissipation 

In this section we study the evolution of a planetary system which starts from a non-resonant 
orbit and evolves under the influence of the dissipation laws ([1]) or ([2]) . 

In Fig. [1^ we present the family of circular orbits and the families of 2/1 and 3/2 resonant 
periodic orbits, in their projection on the X10X20 plane. The family of circular orbits brakes 
at the 2/1 resonance, where a gap appears and we have a bifurcation of two families of 2/1 
resonant elliptic periodic orbits. The same phenomenon is repeated on the circular family at the 
3/2 resonance, where a 3/2 resonant family of elliptic periodic orbits bifurcates (Hadjidemetriou, 
2006). In Fig. [TJd we present a detail of Fig. [T^, where we show the family of non resonant 
circular orbits and the bifurcation of one branch of the 2/1 resonant family. On the circular 
family the 4/1, 3/1 and 5/2 resonant circular orbits are indicated. These orbits are for the 
masses rriQ = 0.9978, mi = 0.0008 and m2 = 0.0014. The family of circular orbits in Fig. [TJj 
is stable except close to the 3/1 resonance. The 2/1 resonant family in this figure is stable, for 
mi < m,2, but a small unstable region appears when m,i > m,2- These will be discussed in the 
following. 
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Figure 1: (a) The family of circular (non resonant) periodic orbits and the families of 2/1 and 
3/2 resonant periodic orbits close to the 2/1 gap. (b) The non resonant circular family and the 
continuation to the 2/1 resonant family after the 2/1 gap. The circular orbits at the 4/1, 3/1 
and 5/2 resonances are indicated. 



2.1 The dissipation law R 
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We start our study by considering a non resonant planetary system with orbital eccentricities 
ei = eio, 62 = 620 and a fixed value of the ratio of the planetary frequencies ni/n2- The initial 
values for the phases M (mean motion) and co (longitude of pericenter) for both planets are 
taken equal to 0° or 180*^, but it was verified that the evolution is similar for any other initial 
values. 



2.1.1 Case I: Starting from ni/n2 ~ 2.8 

We consider a planetary system with initial eccentricities eio = 0.20, 620 = 0.20 and aio = 1.0, 
020 = 2.0, corresponding to ni/n2 ~ 2.8, and for the angles we have taken the initial values 
Ml = 0, M2 = 0, wi = 0, W2 = 180*^. This is a non resonant orbit, whose position could be 
represented in the Fig. [It by a point far from the circular family. The evolution of this system 
in the eccentricity space, under the dissipation law ([TJ with z^ = 5, is presented in Fig. [2j In 
panel (a) we present the first stage of the evolution, where the system is attracted to the family 
of circular orbits, keeping the ratio ni/n2 close to its initial value (ni/n2 ~ 2.8). This motion 
is quite irregular. In panel (b) we present the second stage of the evolution, namely, after the 
trapping of the system in a periodic orbit of the circular family (e^ ~ 0). The system follows a 
route along the circular family, which is also presented, and is attracted to an orbit close to the 
2/1 resonance (third stage), just before the 2/1 gap (see Fig. [TJa). 

In Fig. [3] we present the evolution of the system in time. We use use a logarithmic scale for the 
time to show better the three stages of evolution. During the first stage the semimajor axes, 
and, subsequently, the ratio ni/n2 are almost constant. Up to this point, there is an agreement 
with the results of a first order analysis given by Beauge et al (2006). During the second stage 
the variation of semimajor axes causes a decrease of the ratio ni/n2 up to the resonance 2/1. 
From that point on the ratio ni/n2 remains almost constant and the semimajor axes increase 
(third stage). The eccentricities (panel (b)) decrease at first, while the system moves towards 
the circular family . The eccentricities are finally trapped close to zero (second and third stage 
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Figure 2: The evolution of the system (in the eccentricity space) under the dissipation law ([T]), 
starting from ni/n2 = 2.8 and ei = 0.20, 62 = 0.20. (a) The system is attracted to an orbit on 
the circular family, (b) After it reaches the circular family, the system moves smoothly on the 
circular family until it reaches an orbit close to the 2/1 resonance and is trapped there. 
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Figure 3: The evolution of the planetary orbital elements associated to the system evolution 
presented in Fig. [2J The distinction between the three stages of Fig. O) is clearly seen, (a) 
Evolution of the semimajor axes, (b) Evolution of the eccentricities, (c) Evolution of the ratio 
ni/n2. 
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Figure 4: (a) The first stage of the evolution of the system (in the eccentricity space) under the 
dissipation law R = — 10~^('U — Vc), starting from nxjn^ ~ 3.5 and ei = 0.20, 62 = 0.20. The 
system is attracted to an orbit on the circular family, (b) The second stage of the evolution : the 
system moves smoothly on the circular family, (c) The system is trapped in a chaotic attractor 
close to the circular orbit at the 3/1 resonance. 

of evolution) . 



2.1.2 Case II: ni/n2 ~ 3.5 

We consider a planetary system with initial conditions eio = 0.20, 620 = 0.20, aio = 1.0, 
020 = 2.3, corresponding to ni/n2 ~ 3.5, and angles Mi = 0, M2 = 0, wi = 0, UJ2 = 180*^. 
We study the evolution of this system in the same way as in the previous section, under the 
Stoke's dissipation R = — 10^^(iT — Vc)- The evolution of this system in the eccentricity space 
is presented in Fig. [H In panel (a) we present the first stage and in panel (b) the second stage 
of the evolution, similarly to the previous case. But now, the system is finally attracted to a 
chaotic attractor close to the 3/1 resonance, possibly associated with the unstable region at the 
3/1 resonance (panel (c)). The evolution of the semimajor axes, the eccentricities and the mean 
motion ratio is presented in Fig. [5] and shows similar behaviour as in the previous case. We 
remark the irregular oscillations of the eccentricities near to zero after the system is attracted 
in the chaotic attractor. 

In order to check the universality of the previous results, we repeated the study presented in Fig. 
m for different initial conditions in the eccentricities, keeping the initial ratio nxjn^ ~ 3.5. The 
results are presented in Fig. [6^: In the whole region in the eccentricity space < Cj < 0.30, the 
system is attracted first to the same periodic orbit on the circular family (at ni/n2 ~ 3.5), and 
after this the evolution is the same as that described in Fig. HI We see that this region in the 
eccentricity space is the hasin of attraction towards the chaotic attractor at the 3/1 resonance 
of Fig. Ht. 

We also repeated the study for different values of the planetary masses and verified that the 
evolution is qualitatively the same, as that described above. In figure[6)3 we present the evolution 
of the mean motion ratio for different mass ratios (but for the same total planetary mass). Indeed 
we obtain the same feature of the decreasing ni/n2 but the rate of decreasing, and, subsequently, 
the moment of the final capture in the resonance, depends on the mass ratio. Instead, if we 
change both planetary masses, preserving the mass ratio, the decreasing rate of ni/n2 is not 
affected (see Fig. [6t). 
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Figure 5: The evolution of the elements of the planetary orbits presented in figure [H The 
distinction between the three stages of Fig. HJ) is clearly seen, (a) Evolution of the semimajor 
axes. The semimajor axes increase steadily after the system is trapped close to the 3/1 resonance, 
(b) The evolution of the eccentricities. They remain close to zero even after the trapping to the 
chaotic attractor. c) The evolution of the ratio ni/n2- 
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Figure 6: (a) The evolution of the system, as in Fig. HI for different initial eccentricities, (b) 
The evolution of ni/n2 for different planetary mass ratio values and mi +m2 = 0.0015. (c) The 
evolution of ni/n2 for different mass values and m2/mi = 1.75. 
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Figure 7: The evolution of the system (in the eccentricity space) under the dissipation law 
R = —10~^{v — v'c/r^''^), starting from ni/n2 = 3.5. (a) Starting eccentricities: eio = 0.01 
620 = 0.01. (b) The same evolution as in panel (a) and, additionally, for cio = 0.10 620 = 0.10 
and eio = 0.20, 620 = 0.20 . In all cases the system is attracted to the same asymmetric periodic 
orbit. 



2.2 The dissipation law R 
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In this paragraph we study the evolution of a system under the dissipation law ([2]). We work 
as in the previous section, starting from a planetary system with elements 0.01 < Cjo < 0.2, 
i = 1,2, and aio = 1.0, 020 = 2.3, corresponding to ni/n2 ~ 3.5. For the angles we have taken 
the initial values Mi = 0, M2 = 0, wi = 0, uj2 = ISO'^, but, as it has been mentioned, the 
evolution is similar for any other values of the angles. A presentation of the evolution is given 
in figures [7| and El 

In panel (a) of Fig. [7] we present the evolution for eio = 0.01 and 620 = 0.01. This means that 
the system is a non resonant system with ni/n2 ~ 3.5, with very small eccentricities, almost 
on the family of circular orbits mentioned before and is close to the 3/1 resonant orbit on this 
family. We note that a small unstable region appears on the circular family, at the 3/1 resonance, 
and from the critical points at the ends of this unstable region there bifurcate four symmetric 
families of 3/1 resonant elliptic periodic orbits, one stable and the rest unstable (Voyatzis and 
Hadjidemetriou 2006; Voyatzis, 2008). The ratio ni/n2 decreases rapidly up to the 3/1 resonance 
and from that point on the system follows the path along the stable symmetric family of 3/1 
resonant periodic orbits. However, this family becomes unstable at a certain point, and another 
3/1 resonant elliptic family of asymmetric periodic orbits bifurcates from this critical point. The 
system changes route and follows the stable asymmetric family, until it is trapped at a periodic 
orbit on this family, with relatively large eccentricities. In panel (b) of Fig. [7] we present the 
same evolution starting from different initial eccentricities. All orbits are attracted to the same 
periodic orbit as above. In Fig. [8] we present the evolution of the system in the elements of 
the orbit. In panel (a) we show the evolution of the semimajor axes. They decrease at first, 
and then they remain constant. In panel (b) it is shown that the eccentricities increase at first 
and then they are stabilized to the constant values ei = 0.419, 62 = 0.427. Consequently, the 
evolution is attracted to a particular periodic orbit which is asymmetric as it is clearly indicated 
by the evolution of the angle Alo of the line of apsides of the two planetary orbits. Initially 
the angle Alo oscillates around the value 180° that corresponds to the motion on the symmetric 
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Figure 8: The evolution of Fig. [7^, in the elements of the orbit, (a) The evolution of the 
semimajor axes, (b) The evolution of the eccentricities, (c) The evolution of the planetary 
apsidal difference Acj starting close (Acjo = 0.01) and far (Aejo = 0.1) from the periodic orbit; 
the transition from the symmetric to the asymmetric part of the 3/1 family is clearly seen. 
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Figure 9: The asymmetric periodic orbit to which the system is trapped in the evolution of Fig. 
[71 in (a) the inertial and (b) the rotating frame. 



family and finally tends to the value Aw = 87°. 4, which indicates clearly the asymmetry of the 
orbit. 

This asymmetric orbit, which is the attractor of the evolution, is shown in Fig. [H In order 
to stop the planetary migration process and obtain an attractor orbit, Lee and Peale (2002) 
artificially introduce to the system a sufficient amount of eccentricity damping. In our case the 
migration stops naturally if the dissipation law ([2]) is taken into account. 

Finally, we present in Fig. [TOb the evolution of the ratio ni/n2, under the dissipation law 
R = —W~'^{v — v'c/r^''^), for u = 5 and u = Q. We found that in almost all cases the system is 
trapped at the 3/1 resonance, as shown in Fig. [71 We found only one special case, corresponding 
to z/ = 5 and eio < 0.01, 620 < 0.01, where the system is not trapped to the 3/1 resonance, as 
the ratio ni/n2 decreases, but goes up to the 2/1 gap. From that point on, it is trapped and 
moves on the stable family of 2/1 resonant elliptic periodic orbits, until finally it is attracted to 
a 2/1 resonant periodic orbit on the 2/1 family, with large eccentricities. This is shown in Fig. 
[TOb. 
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Figure 10: (a) The evolution of Fig. [3 in the ratio ni/n2 starting from ni/n2 = 3.5. For z^ = 6 
the system is trapped in the 3/1 resonance but for z^ = 5 the system is trapped in the 2/1 
resonance, (b) The trapping in the 2/1 resonance of panel (a) in the eccentricity plane. 

3 The evolution close to resonant orbits under dissipation: 
2/1 resonance 



3.1 The dissipation law R 
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In this section we study the evolution of a planetary system under the dissipation law ([T]). We 
consider orbits which start in a region close to the families of 2/1 resonant symmetric periodic 
orbits. Since the structure of the families of periodic orbits is essentially different between the 
cases mi < m2 and ttt-i > ?tt.2, we study both cases separately. A detailed study of the families 
of periodic orbits at the 2/1 resonance is given in Voyatzis and Hadjidemetriou (2005), Beauge 
et al (2006), Voyatzis et al (2009). We mention that similar dynamics with that presented below 
is obtained for the 3/1 resonance too. 



3.1.1 The inner planet has a smaUer mass than the outer planet. 

We used the mass values mi = 0.0008 and m,2 = 0.0014. We start in a region close to the 2/1 
resonant family presented in Fig. [T^. There are two resonant symmetric families, separated by 
the 2/1 gap. In Fig. [TTb we present these two families in the eccentricity space 6162, denoted by 
I and II. In family I both planets are at perihelion (Acj = 0) and in family II the inner planet 
is at perihelion and the outer at aphelion (Acj = I8O''). The 2/1 gap of Fig. [1^ appears in 
this diagram as a small gap at ei ?« 0, 62 ~ 0. The family I is all stable. The family II starts 
as unstable, for small eccentricities, until a gap appears, corresponding to close collision orbits 
between the planets. After this collision area, the family II is stable up to large eccentricities. 

We study the evolution of the system by starting from the two points indicated in Fig. [TTk . The 
first point is on the family, with high eccentricities (ei = 0.80, 62 = 0.40). The second point is 
obtained by shifting the value ei by 0.14 and keeping the value of 62 the same. We indicate the 
two evolution cases as A and B. The evolution, for z/ = 7 is shown in Fig. Illb . In the case A, 
the system moves along this family, with decreasing eccentricities (in absolute values), until it 
is trapped to a periodic orbit with almost zero eccentricities, close to the 2/1 resonance. In the 
evolution of orbit B, the system also follows a path along the family, which acts as a guiding 
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Figure 11: (a) The two families of 2/1 resonant periodic orbits and the two starting points for 
the evolution of the system under dissipation, (b) The evolution of the elements of the planetary 
orbits under the dissipation law R = —10~'^{v — Vc), starting on the family I of the 2/1 resonant 
periodic orbits (case A) and far from it (case B). 

line, with oscillations around the family that decrease in amplitude and tends to the same zero 
eccentricity orbit as in the case A. 

In Fig. [12] we present the evolution of Fig. Illb . in the elements of the orbit. In panel (a) we 
show the evolution of the semimajor axes of the two planets. They decrease at first and then are 
stabilized, as the system is trapped at the 2/1 resonant periodic orbit with zero eccentricities. 
In panel (b) we show the evolution of the eccentricities. They decrease (in absolute values) at 
first and stay at small values after the trapping to the zero eccentricity orbit mentioned above. 
The ratio ni/n2 stays close to the 2/1 value (panel (c)). 

In the previous study, we considered a system close to the family I of 2/1 resonant periodic 
orbits of Fig. [TTk . We study now a system that starts close to the family II in this figure. 
The evolution is shown in Fig. [131 The starting point is close to ei = 0.45, 62 = 0.6 on the 
family II. The system moves along the family, with decreasing eccentricities, until it is trapped 
temporarily at a point on the family, with smaller eccentricities, close to the collision area. The 
detail of this trapping is shown in the separate plot of Fig. [T3j However, after this temporary 
capture, the system is disrupted due to the existence of strong chaos nearby the collision region. 



3.1.2 The inner planet has larger mass than the outer planet 

The fact that the inner planet is more massive than the outer planet results to an important 
qualitative diff'erence for the family I of resonant periodic orbits: A small unstable region appears 
when mi > m2, which increases as the ratio mi/m2 increases (see Beauge et al, 2006; Voyatzis 
et al, 2009). For both critical points formed, there bifurcate two asymmetric families of 2/1 
resonant periodic orbits In the following we use the masses mi = 0.0015, m,2 = 0.0010. For 
these values the two bifurcating asymmetric families meet and form a single family. 

In Fig. [T^ we study the evolution starting from a point on the symmetric family, with high 
eccentricities. The system moves along the stable symmetric family until it meets the unstable 
region, where it changes route and follows the asymmetric family until this latter family meets 
again the symmetric family. From that point on it follows the path of the stable symmetric 



1.6 -\ 



(B) 



(A) 



(B) 



(A) 



2x10' 4x10' 

time (t.u.) 



(a) 



6x10' 



ii 



f^::^ 



/ 



/ 



Case (B) 



2x10' 4x10' 

time (t.u.) 



(b) 



6x10' 







(c) 


2.06- 




^ 




(^' Jl^ 


r 


2.04- 


/ 




2.02 - 


Wiir "■ 






1 ' 





0x10° 



2x10 4x10 

time (t.u.) 



Figure 12: The evolution of Fig. Illb . in the elements of the orbit. Both the evolution along the 
family (case A) and close to the family (case B) are shown, (a) The evolution of the semimajor 
axes, (b) The evolution of the eccentricities (only case B is shown), (c) The evolution of n\lni. 
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Figure 13: The evolution of the system on the family II of the 2/1 resonant periodic orbits 
(see Fig. [TTb). The starting point is at high eccentricities. The eccentricities decrease and the 
system is trapped in a bounded chaotic region close to the collision area, before widespread 
chaos appears. A detail of the evolution close to the temporary chaotic trapping is shown in the 
separate plot. 
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Figure 14: (a) The evolution along the 2:1 resonant family of periodic orbits for rfi\ = 0.0015, 
7712 = 0.0010, starting from high eccentricities. The evolution follows the stable parts of the 
family, (b) The evolution of Auj, which indicates the passage of the system through different 
symmetric and asymmetric configurations. 

family and is finally trapped at a 2/1 resonant periodic orbit with zero eccentricities. The 
evolution is qualitatively similar for u = 5 and z^ = 6 in the dissipation law ([T]). In Fig. 114b we 
present the evolution of the angle Aw between the line of apsides. The transition from symmetry 
Auo = to asymmetry (AtJ ^ 0^, 180°) and again to symmetry, AtJ = and then to Au: = 180°, 
is clearly seen. 



3.2 The dissipation law R = -10^^(7;- Vc/r^^'^) 

In this section we consider the dissipation law ([2]) and repeat the same study as in section 13. H 
starting from two orbits on the family of 2/1 resonant periodic orbits, one with large eccentricities 
and the other with small eccentricities. We consider two cases: (1) mi < m2, using the masses 
mi = 0.0008, m,2 = 0.0014 and (2) m,i > m,2, using the masses mi = 0.0015, m,2 = 0.0010. In 
case 1, the family of periodic orbits is given in Fig. [TTk . family I, and in case 2 the family of 
periodic orbits is given in Fig. [T^ . 

The evolution is very different from that of Fig. [TTJ Starting from large eccentricities, the 
evolution is towards smaller eccentricities and starting from small eccentricities, the evolution 
is towards larger eccentricities. In both cases, the system is trapped to a 2/1 resonant periodic 
orbit on the corresponding family. The results are presented in Fig. [TSl In panel (a) we show 
the evolution when we start close to the family I of Fig. [TTk . We start both with large and with 
small eccentricities (absolute values), and in each case we consider two starting points, one on 
the resonant family and one close to it. The system is attracted to the same orbit on this family, 
with ei = 0,442, 62 = 0.137. If we start on the family, the system moves always on the family. 
If we start close to the family, then the system oscillates around the resonant family of periodic 
orbits, with decreasing amplitude, and is finally trapped to the periodic orbit mentioned above. 
We remark that the attracting orbit of Fig. [T5b is the same as the attracting orbit of Fig. [TOl 
In panel (b) we do the same work, starting close to the family of Fig. [T^ . In this latter case, it 
is m.1 > 7712 and the stable family includes an asymmetric branch. The evolution is the same as 
before, moving always along the stable family. In the present case the orbit to which the system 
is trapped is an orbit on the asymmetric branch of the 2/1 resonant periodic orbits. 
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Figure 15: The evolution of the system (in the eccentricity space) under the dissipation law 
R = —lQ~^{v — Vc/f^ ): along a 2/1 resonant family, starting on the family and close to the 
family, both from small and from large eccentricities (absolute values). The system is trapped 
to a 2/1 resonant periodic orbit, (a) The evolution close to the family I for mi < m2 (see Fig. 
[TTb). (b) Similar evolution as in panel (a), close to the family I for mi > m2 (see Fig. [Hk). 

4 The evolution close to 3/2 resonant orbits under dissipation 

In this section we present the evolution close to the 3/2 resonance, under three different dissi- 
pation laws. The results are presented in Fig. [TGk.b.c. In all panels the family of 3/2 resonant 
symmetric periodic orbits is presented, in the e\ 62 space (the same family in the plane xi X2 is 
presented in Fig. [1^). Along the family, the stability changes and after a certain point insta- 
bility appears. At the critical point, a family of 3/2 resonant stable asymmetric periodic orbits 
bifurcates (Michtchenko et al, 2006). In all cases we used the planetary masses mi = 0.0008, 
m2 = 0.0014. 

In panel (a) we study the evolution under the law ([1]), particularly R = — 10~^('U — Vc). The 
starting point is at large eccentricities and the system evolves to zero eccentricities. This evo- 
lution is similar to that obtained for the 2/1 and 3/1 resonance but now the semimajor axes 
decrease slightly at first and then they are stabilized to fixed values and the system is trapped 
in a particular circular orbit. 

In panel (b) we present the evolution of the system under the law (l2|), particularly R = — 10~^(u— 
Vc/r^''^)- We start from two points, one with large and one with small eccentricities, as in the 
case of Fig. [T5k . The evolution is similar as in this latter figure: We find that the system 
is attracted to a periodic orbit with ei = 0.187, 62 = 0.196, on the 3/2 resonant family, by 
decreasing or increasing eccentricities, respectively. 

In panel (c) we study the evolution under the law R = — 10~^(i; — u^/r^'^), which differs from 
the law used in panel (b). In this case the velocity of rotation of the protoplanetary nebula is 
stronger than the Keplerian circular velocity Vc close to the sun and decreases more rapidly for 
larger distances, compared to the dissipation law of panel (b). We start from small eccentricities 
and we note that the eccentricities increase, until the critical point where instability starts. From 
that point on, it follows the stable 3/2 asymmetric family, up to a certain point before chaos 
appears and the system is disrupted. 
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Figure 16: The evolution close to the 3/2 resonance, under three different dissipation laws, (a) 
R = _io-6(i;_^;-;). (b) R = -10-6(^-t;-;/rV2) (c) R = -10-^(u-?;Vr3/2). The starting point 

and the direction of the evolution along the family of periodic orbits are indicated in all cases. 

5 Discussion 



In this work we studied the evolution of a planetary system with two planets moving in the 
same plane, under the action of dissipative forces. The main study was made for the Stokes 
like force given by Eq. ([T|), for different values of v. This force implies that the two planets 
move inside a protoplanetary nebula that rotates with a differential Keplerian circular velocity 
Vc- We also considered a variation of this dissipative force given by eq. ([2]) and, for the 3/2 
resonance, also the force R = — 10~^(i/ — Vc/f'^ )■ In these last cases the rotational velocity 
of the protoplanetary nebula deviates from the circular Keplerian velocity, being larger close to 
the sun and tending to zero at large distances. 

An important result that was obtained in the present study is that the families of periodic 
orbits play a crucial role in the evolution of the system: These families, both stable or unstable, 
determine the structure of the phase space. It was shown that the stable families provide the 
routes along which the system evolves, while the unstable families define the regions that the 
system avoids in its evolution. The evolution may he different for different dissipation laws, hut 
in all cases it is the families of periodic orhits that guide the system in its evolution, either to 
smaller or to larger eccentricities. The particular dissipation law determines how the system 
moves along the families and its final trapping to a resonance. 

We studied the evolution of a planetary system considering two different cases. 

Case 1: The planetary system starts far from a periodic motion (either resonant or non resonant). 
Case 2: The planetary system starts close to a resonant periodic motion. In particular, we 
studied the 2/1 and 3/2 resonances. 

The evolution depends on the particular dissipative force. 

Case 1 : 

We considered first the Stokes like law of eq. ([T]). We found that in all cases the system is 
attracted to an orbit on the circular family of periodic orhits, following an irregular route with 
almost the same value of the planetary frequencies nxjui as the starting value. After this, 
the system evolves along the circular family, in a smooth way, to the direction where nxjn^ 
decreases, until it meets a significant resonance, for example 2/1 or 3/1. From that point on, 
the system is trapped in a chaotic attractor or in a resonant orbit, respectively, with very small 
eccentricities (see figures [2l [3|) . The semimajor axes increase slowly and then are stabilized at 



the values corresponding to the resonant attractor. This evolution is the same for a whole region 
in the eccentricity space, which is in fact the basin of attraction to the above mentioned chaotic 
attractor (Fig. [6]) for the 3/1 resonance or to a circular resonant orbit for the 2/1 resonance. 
Thus we get a mechanism of trapping to a resonance with very small eccentricities. 

A different evolution appears for the dissipation law ([2]). The system moves towards a 3/1 
resonant orbit on the circular family and then follows the stable family of elliptic 3/1 resonant 
periodic orbits that bifurcates from the resonant circular orbit. It evolves along this family, 
with increasing eccentricities, and if the family becomes unstable, it changes route and follows 
the stable asymmetric family that bifurcates from this critical point. This process provides a 
mechanism of trapping to an asymmetric 3/1 resonant periodic orbit with large eccentricities. 

Case 2 : 

We start with the Stokes like eq. ([T]), as in case 1. We study first the evolution close to the 
2/1 resonance. In all cases, the evolution follows the stable parts of the family of 2/1 resonant 
periodic orbits, with decreasing eccentricities (see figures [Tib . [T4b for the 2/1 resonance and 
Fig. [T6k for the 3/2 resonance). The semimajor axes are stabilized to fixed values, after a 
slow decrease. If, however, there is a chaotic region along the family, then the system develops 
chaotic motion (unless the dissipation process stops earlier) with decreasing eccentricities (Fig. 
[T3]) . This evolution is the same either for a starting orbit on the 2/1 family or far from it. 

The evolution is different for the dissipative force given by eq. ^. Depending on the starting 
point, the eccentricities decrease or increase and the system is finally trapped to a resonant 
periodic orbit with nonzero eccentricities (see Figll5lfor the 2/1 resonance and Fig. 116b for the 
3/2 resonance). 

Finally, a different evolution appears for the dissipative force R = —10~'^{v — v'c/r^''^). Contrary 
to the previous two cases, the eccentricities in this case increase, following the stable parts of 
the resonant family. Chaotic motion may appear if an unstable region is met on the family (see 
FigJlGb for the 3/2 resonance). 

We see that the final outcome, whether the system is trapped to small or large eccentricities, 
depends critically on the particular dissipation force, and in particular on how much the rota- 
tional velocity of the protoplanetary nebula deviates from the Keplerian circular velocity. In 
almost all cases the system is trapped to a resonant orbit. 
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